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Abstract: We present analytical results to guide numerical simulations with Wilson 
fermions in higher representations of the colour group. The ratio of A parameters, 
the additive renormalization of the fermion mass, and the renormalization of fermion 
bilinears are computed in perturbation theory, including cactus resummation. We 
recall the chiral Lagrangian for the different patterns of symmetry breaking that can 
take place with fermions in higher representations, and discuss the possibility of an 
Aoki phase as the fermion mass is reduced at finite lattice spacing. 
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1. Introduction 

The nonperturbative dynamics of asymptotically free gauge theories with matter 
fields transforming according to higher dimensional representations of the underlying 
gauge group is a topic of current research interest. Physics beyond the Standard 
Model (BSM) offers a new arena for the use of these theories, which are expected to 
develop a nonperturbative infrared fixed point for a very low number of flavors Jl|, 0] . 
Because they minimize the tension with the electroweak precision constraints, some of 
these theories are excellent candidates for the dynamical breaking of the electroweak 
symmetry of walking technicolor type, which were first introduced in Refs. |3], [5], 
H, 0, §] . Since the number of flavors needed to get near the conformal fixed point is 
small, the associated models have been termed minimal walking technicolor |], ||, §]. 
By walking one refers to the fact that the running coupling decreases much more 
slowly with the reference energy scale than in the case of QCD-like theories. Yet, 
another interesting physical application of the study of the phase diagram of strongly 
coupled theories is to provide the theoretical landscape underlying the unparticle 
physics world |T0, TT]. The theory landscape was provided in Ref. JT^] where it was 
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shown that the fraction of asymptotically free gauge theories developing an infrared 
fixed point is quite large. Studying their phase diagram is a fundamental step if 
these theories aspire to become realistic candidates for BSM physics. Insight into 
the phase diagram of such theories has been recently provided by a proposed all-order 
beta function for any number of colours and for any representation [ I3| . Moreover, in 
the limit of a large number of colours, the planar orientifold planar equivalence relates 
theories with fermions in higher representations to supersymmetric theories It 
provides interesting predictions that deserve nonperturbative investigations fl5], |16] . 
The necessity to study the large-iV limit makes these theories more expensive to 
study numerically. Finally, understanding the strong dynamics that governs the 
low-energy behaviour of such theories is an interesting problem per se. 

Lattice is a privileged tool for exploring the nonperturbative dynamics of strongly 
interacting theories, but Monte Carlo simulations of these theories can capture the 
interesting dynamical features only if the full fermion determinant is taken into ac- 
count in the Boltzmann weight used for generating gauge configurations. So far 
only limited experience has been gathered from numerical simulations with dynam- 
ical fermions beyond QCD ||T7| , 18] . In the light of recent algorithmic progresses in 



simulating quantum field theories with dynamical fermions [19, 20, 21, 22, 23, 24 



numerical studies with fermions in higher-dimensional representations are now a re- 
alistic target (see Ref. |1| for early work in this direction). As a preliminary work, 
which provides guidance for large-scale simulations, we present here an investigation 
of the space of bare couplings by analytical tools, such as perturbation theory and 
chiral Lagrangians. 

Perturbative results are useful to understand the behaviour of the lattice theory 
as the continuum limit is approached: on one hand they provide a connection between 
the lattice results and their continuum counterparts; on the other hand they offer 
some quantitative support in choosing the bare parameters in the early stages of 
numerical simulations. Precision studies in QCD have shown sizable discrepancies 
between perturbative and nonperturbative computations at the values of the bare 
parameters that are currently accessible. Assessing the accuracy of perturbation 
theory for theories with fermions in higher-dimensional representations is beyond 
the scope of this paper and will be deferred to future publications. Instead we 
shall supplement perturbative calculations with sensible assumptions, that we discuss 
below, in order to dictate the choice of the values of the bare parameters for first 
numerical investigations. 

In this work perturbation theory is used to determine the dependence of the lat- 
tice spacing on the bare lattice coupling, the ratio of lambda parameters Aj^g/Ai at , 
the critical value of the bare mass m c and critical hopping parameter k c , and the 
renormalization constants for fermionic bilinears. The gauge action considered is 
the simple plaquette action, and the fermion action is the unimproved Wilson ac- 
tion. This simple choice provides a concrete example for performing the perturba- 



- 2 - 



tive calculations, and matches the existing, and forthcoming numerical simulations. 
Four specific examples of lattice theories with Wilson fermions are compared below. 
Results for quenched QCD (QCDO), and for QCD with two flavours of dynamical 
fermions (QCD2) are known both in perturbation theory, and from non-perturbative 
computations. They are briefly summarized in this work in order to set the frame- 
work for our computations, and to assess the accuracy of perturbation theory. We 
then present the generalization of the perturbative calculations to arbitrary repre- 
sentations, and analyze in detail their implications for two theories that are good 
candidates for BSM phenomenology, namely the SU(3) gauge theory with rif = 2 
flavours in the two-index symmetric representation (Tl), and the SU(2) gauge theory 
with two flavours in the two-index symmetric representation (T2). 

The chiral Lagrangian describing the dynamics of the light Goldstone bosons is 
analyzed in order to clarify the structure of the phase diagram that is likely to be 
revealed by numerical simulations for small quark masses. We discuss in particular 
the theories with fermions in higher representations introduced above. Note that the 
same approach can be readily applied to other theories, like e.g. gauge theories with 
fermions in the two-index antisymmetric representations, that are interesting for 
numerical tests of the planar orientifold equivalence. Even though we do not discuss 
these theories explicitly here, our conclusions can be specialized in a straightforward 
manner. 

The paper is organized as follows. In Sect. we discuss the perturbative results 
describing asymptotic scaling, we compute the ratio of A parameters, and discuss 
the approach of the continuum limit. Sect. |3] reports some useful results at one- 
loop in perturbation theory. We first consider the renormalization of the bare mass 
for Wilson fermions; the critical value of the hopping parameter k c is computed 
both up to two loops, and using the so-called cactus dressing to resum a particular 
class of tadpole diagrams [25]. Similarly we present results for the renormalization 
constants for fermion bilinears. Finally in Sect. |] we discuss the form of the chiral 
Lagrangian that describes the low-energy dynamics in theories with fermions in 
arbitrary representation and the possible phase structure as the quark mass is lowered 
at finite lattice spacing. 

Numerical simulations of the theories with fermions in higher representations are 
deferred to further publications. 



2. Scaling 

2.1 Ratio of A parameters 

The P function encodes the dependence of the lattice spacing a on the bare coupling 
constant g . In mass-independent renormalization schemes, the lattice spacing is 
uniquely determined by the bare coupling, according to the renormalization group 
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equation: 



AatG?o) = -a^p , (2.1) 

0(1 g R fixed 

where the subscript lat refers to the lattice scheme which is being considered here. 

For a generic gauge theory, with gauge group SU(iV) and rif fermions in a given 
representation R of the colour group, the two-loop computation in perturbation 
theory yields the familiar expression: 
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where T R yields the normalization of the generators, and C 2 (R) is the quadratic 
Casimir, both in the representation R. Factors of C 2 (A), the quadratic Casimir in 
the adjoint representation, arise because of gluon loops, and do not change as the 
fermionic representation is varied. Note that the first two coefficients of the (3 function 
are universal and depend neither on the regularization nor on the renormalization 
scheme. For N = 3 and fermions in the fundamental representation, the expressions 
above reduce to the usual values of Po, Pi- Tables for the group-theoretical factors 
for the representations considered in this work are reported in App. |A]. 

The asymptotic behaviour of a(go) is obtained by integrating Eq. (|2.1|), and the 
scale Ai a t is the integration constant that appears in this procedure. Following the 
notation in Ref. we write: 
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The ratio of A-parameters defined in different renormalization schemes is ob- 



tained from the one-loop relation between the coupling constants f27| , p28| , g9[ Pq| . 
In particular the running coupling in the MS scheme is related to the lattice bare 
coupling via: 



-1/2 



9o, 



;2.6) 



1=1 



where Ao is the bare gauge-fixing parameter in the lattice formulation. This relation 
can be obtained e.g. using the background field technique 0, [3"2]1 . The first 
coefficient (/ia, Ao) for matter fields in the fundamental representation can be 



found in the literature, see e.g. Ref. [29, 26 1. It has the generic form: 



2 {1) {^a, Ao) 
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Representation 










SU(2), n f 


= 




0.0464389 


0.00181793 


-0.277412 


SU(3), n f 


= 




0.0696583 


0.00409035 


-0.468201 


SU(2), n, 


= 2, 


fund 


0.0612149 


0.00307445 


-0.454809 


SU(3), n/ 


= 2, 


fund 


0.0612149 


0.00307445 


-0.454809 


SU(3), n f 


= 2, 


2S 


0.0274412 


-0.00259323 


-0.401241 


SU(2), 


= 2, 


2S 


0.0126651 


-0.00160406 


-0.223844 



Table 1: Perturbative coefficients appearing in one-loop perturbative computations. 



and the ratio of A-parameters is obtained from the coefficients in Eq. ( |2.7| ) as: 

A lat /A M g = exp[/ /(2/? )]. (2.8) 

Having already written the coefficient f3$ for a generic representation in Eq. ( |2.3| ), 
the expression for the finite part of the one-loop contribution in Eq. ( j2.6[ ), Iq, is 
the only ingredient needed in order to convert the A-parameter. The coefficient Iq 
is obtained by inspecting the one-loop diagrams that contribute to Eq. (|2.7|) . The 
group-theoretical factors need to be changed in order to take into account the new 
fermionic representation, while the numerical factors that arise from the integration 
over the lattice momenta remain unchanged. For a generic representation R, we 
obtain: 

*o = ttAv h 2 ^ 2 C 2 (F) - 3.54958342046 C 2 (A) + 

+ 1.057389936 T R n f ] , (2.9) 

where the only dependence on the fermionic representation is encoded in the last 
term in the sum on the RHS of Eq. fl2.9| ). Explicit values for some representations of 
interest are summarized in Tab. |I|. The well-known values for the quenched SU(2) 
and SU(3) theories, and for QCD with two flavours of fundamental fermions are 
reported in order to show explicitly the differences in the perturbative coefficients 
as we introduce matter in higher representations. As a non-trivial check, we can 
specialize Eq. ( |2.9|) to the case of M = 1 SYM, which corresponds to n/ = 1/2 
flavour of fermions in the adjoint representation. Using the group-theoretical factors 
reported in Tab. |2| in Appendix |A], our expression reproduces Eq. (24) in Ref. p^ . 

The ratios of A-parameters are easily obtained from the values in Tab. [l|. The 
coefficients in the first lines of the table reproduce the known results: 

19.82, (2.10) 
28.81, (2.11) 
41.05, (2.12) 



7 *-MS/^'lat|su(2),n / =0 
^■Ms/-^lat|su(3),n/=0 
A-MsMlat I SU(3),n / =2,fund 
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see e.g. Refs. pq , 29, 26]. The last two lines yield the new results for the higher 



representations that we want to consider in this work: 

AMs/Aiatlsu(3),n /= 2, 2 5 = 14 69-59 (2.13) 
AMs/Aiat|su(2),n /= 2,2S = 6884.36. (2.14) 

Values for other representations can be easily deduced from the formulae above. 
Including two flavours of fermions in higher representations induces large variations 
in the ratios of A parameters. This is at odds with the results for fermions in the 
fundamental representation, where adding the effect of fermion loops yields a much 
smaller variation of the ratio. However the large values obtained for fermions in 
two-index representations can be understood by rewriting the ratio /3o/2/o in a way 
which makes the 1/N scaling explicit: 

l 1 - 0.0787969T fl ^f - 0.735484^ 

~ -3.65978 A nf — ■ 2.15 

2A) l ~Ti T RW 



In Eq. ([2.1 5|) we can easily recognize the contributions 0(rif/N) from fermion loops, 
and the contributions from the non-planar gluonic diagrams of 0(1/N 2 ). For fer- 
mions in the fundamental representation Tr = 1/2, and therefore the fermion de- 
terminant yields corrections that are suppressed by n,f/N. For fermions in 2-index 
representations, the normalization of the generators is such that Tr ~ O(N), and 
hence the contribution from the fermion determinant is of the same order of the 
gluon contribution, both in the numerator and the denominator, and therefore large 
variations are found with respect to the pure gauge theory. 



2.2 Perturbative and nonperturbative scaling 

The perturbative results obtained in the previous subsection can be used to sketch the 
scaling of the lattice spacing for theories in higher representations, being well aware 
of the limitations of perturbation theory. The nonperturbative scaling of the lattice 
spacing has been carefully studied for QCD, both in the quenched approximation, 
and for the theory with dynamical quarks in the fundamental representation. For the 
latter theories, the accuracy of perturbative estimates can be assessed by comparing 
numerical and analytical results. As we shall see below, perturbation theory in QCD 
does not yield an accurate description of the scaling of physical quantities. Therefore, 
any result obtained in this framework is bound to be approximate and should be used 
mostly as a guide for forthcoming numerical simulations. 

In order to relate more easily to the notation used in numerical simulations, let 
us introduce the lattice coupling (3 = 2N/gl. We will henceforth use (3 to indicate 
the bare lattice coupling, unless explicitly stated. The asymptotic scaling formula 
reported in the previous subsection yields the value of the lattice spacing in physical 
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units as a function of (3: 



^MS 



A 



lat 



^■MS 



exp 







4iV/3 



(2N(3 /(3) M2 ^ ] [1 + 0(1/(3)] 



(2.16) 



Having computed the ratio of A parameters, the only input that is required is the 
value of Ajyfg. 

For the SU(3) pure gauge theory, Fig. [1] displays the prediction for the lattice 
spacing a in physical units [fm] , computed from two-loop perturbation theory using 
the input from Ref. [[H. [35|]; the curve is compared to the interpolation of the nonper- 



turbative data presented in |36|]. The error band in the figure is simply the error that 
is obtained from propagating the error in the determination of A^ to the value of 
a{(3). As shown by the plot, the perturbative prediction in bare perturbation theory 
underestimates the actual lattice spacing by 30-50% at the values of (3 between 5.8 
and 6.2, where most simulations have been performed so far. The two computations 
agree at large values of /?, as expected when the continuum limit is approached. 



0.14 




Figure 1: Comparison of the lattice scaling in physical units predicted from perturbation 
theory with the nonperturbative results obtained from numerical simulations. The theory 
is pure gauge SU(3). The blue (respectively red) curve represents the perturbative (resp. 
nonperturbative) estimate of the lattice spacing in fm as a function of the lattice bare 
coupling (3. The error on the perturbative estimate comes from the error in the determi- 
nation of Ajy[g. The red curve is an interpolation of the nonpertubative determination of 
the lattice spacing. 

For the SU(3) theory with two flavours of Wilson fermions, the perturbative 
prediction is obtained using the A parameter computed in Ref. [B7], B3|; it can be 
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compared to the value of the lattice spacing recently computed for two values of 
(3 in Ref. f39fl . Again, in the range that is accessible to current simulations, the 
perturbative estimate is smaller by approximately 30%-40%. Such large deviations 
between lattice bare perturbation theory and nonperturbative results are very well- 
known to lattice QCD practitioners P0| . They are reported here in order to have a 
concise summary of the results in QCD, before moving into new territories. 

Figure shows the perturbative estimate for the value of the dimensionless quan- 
tity a((3)Ayfir for QCD with n/ = 0,2 in the same range of /?; for the values of (3 
used in current simulations, one can see that l/(aA^jg) « 20. Given that the hadron 
masses in QCD turn out to be of the order of A^jg, the value of l/(aAj^g) is such that 
lattice artefacts are small, while sufficiently large physical volumes can be reached 
on lattices that have 20-30 points in each direction. Lattice simulations for theories 
beyond QCD need to identify a similar regime in order to avoid large lattice artefacts 
and/or large finite volume effects. 




Figure 2: The perturbative result for the dimensionless quantity a{(3)A.ym as a function 
of (3 for pure gauge SU(3) (blue line), and SU(3) with rif = 2 flavours of Wilson fermions 
(red line). 



Perturbative results can be used for the theories with fermions in higher rep- 
resentations in order to arrive at an educated guess for the value of Ky^a{f3) from 
perturbative scaling, provided a few hypotheses are made in order to identify the 
relevant energy scales. As already mentioned in the introduction, a near-conformal 
behaviour is expected in the theories Tl and T2. The dependence of the coupling 
on the lattice spacing in these theories is characterized by two different regimes. At 
high energies the theories are asymptotically free, and therefore we expect the usual 
logarithmic running of the renormalized coupling. However, as the energy scale is 
decreased, it should reach a value, which we denote A w , where the coupling starts to 
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"walk", i.e. where the coupling is only weakly dependent on the cutoff. The walk- 
ing behaviour should extend for several orders of magnitude, until a lower scale Air 
where the coupling starts running again. Phenomenologically relevant models would 
favour a ratio A W /A IR > 10 3 0. However, it should be noted that the running of 
the coupling constant is scheme dependent, and therefore this ratio should only be 
taken as an indicative value. 

The perturbative values of A^a(P) as a function of f3 are reported in Fig. |3| for 
the theories Tl and T2. If we assume that A^jg ~ Air, and if we further require that 
A w a(/3) ~ 0.1, then numerical simulations of phenomenologically relevant models 
would require Ky^a{f3) < 1CT 5 . The range of (3 in the figure is chosen to yield 
values of Aj^a(P) that saturate the inequality above; if the ratio between the typical 
hadron masses turn out to be large in units of A^g, so that m^^a ~ 1, then higher 
values of (3 would be needed in order to keep both lattice artefacts and finite volume 
effects under control. More quantitative information can only be obtained from 
numerical simulations. As already seen above, in numerical simulations of QCD 




Figure 3: The lattice spacing in unit of A^jg as a function of /3 for the theory Tl (left), 
and T2 (right). 



nonperturbative scaling does deviate from the two-loop predictions by up to 50%. 
We should therefore take these perturbative results with a grain of salt. 

It is worthwhile to emphasize that, while we have expressed everything in units 
of Ajyfg, the absolute value of the scale for these new theories is not known a priori, 
and can only be determined by computing some physical dimensionful quantity. An 
obvious candidate would be the decay constant of the technipion in the chiral limit, 
which is related to the Higgs vev and can be estimated to be F ~ 250 GeV. 

Nonperturbative results that could highlight the near-conformal behaviour can 
be obtained by Monte Carlo renormalization group methods [4^]. However these 
methods require to vary the lattice cutoff over a large interval, while simultaneously 
keeping the finite-size errors under control. If the physical scales of interest are well 
separated, lattice simulations require a very fine resolution, i.e. a large number of 
lattice points, that may not be accessible with present-day computing resources. 
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The finite-volume schemes introduced in Ref. [43]| provide an elegant solution 
to this problem; the renormalization scale is identified with the inverse of the linear 
lattice size, and the evolution of the renormalized coupling is computed in steps, 
changing the scale [i by factors of 2 in each step. The variation of the coupling is 
summarized in the step-scaling function, which yields a precise determination of the 
nonperturbative beta function. Recent results for the theory with nf = 12 fermions 
in the fundamental representation show that this is a promising way to study the 
problem El. 



3. Perturbative renormalization 

In order to perform numerical simulations, preliminary estimates of the critical mass 
and of the renormalization of fermion bilinears are needed. This section summarizes 
some useful computations at one and two loops in perturbation theory for fermions 
in generic representations, that can be used to guide preliminary lattice studies. 
We will consider here the theory defined on the lattice, with Wilson action for the 
fermions, and simple plaquette action for the gauge fields. In the gauge action, the 
link variables are always SU(iV) matrices in the fundamental representation, while in 
the fermionic part of the action, the covariant derivatives are defined through the link 
variables in the generic representation. Feynman rules for perturbative calculations 
are easily generalized. 

Simulating Wilson fermions, the bare mass in the Lagrangian undergoes an ad- 
ditive renormalization, so that the chiral limit is reached for a critical value m c which 
needs to be determined nonperturbatively. Again perturbation theory is useful to 
get some guidance on the initial choice of parameters before embarking in actual 
simulations. 

Following the notation in Ref. [45|1, we write the perturbative expansion of the 



one-particle irreducible two point function at zero momentum as: 

m c (g ) = g 2 ^ + g$P> + . . . (3 .i) 

At one loop, the usual tadpole and sunset diagrams give rise to two contributions, 
Ci and , respectively; for a generic fermionic representation R, these yield: 

= 2C 2 (R) [c[ 1} + 4 1} ] . (3.2) 

At this order in perturbation theory, the additive renormalization of the mass is 
simply proportional to the quadratic Casimir of the fermionic representation, while 
the proportionality constant + c^ 1 is independent of the representation. We can 
therefore use the value in Ref. [ PH : 



"1 ' c 2 



-0.162857058711(2). (3.3) 
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At two loops we need to inspect the structure of the diagrams listed in Fig. 2 of 



Ref . J|5j , and modify the group theoretical factors to take into account the fact that 
the generators appearing in the vertices involving both fermions and gluons are in 
a generic representation, whereas generators inside 3- and 4-gluon vertices are still 
in the fundamental representation. For our present purposes, the result of Ref. [|5 
(Eq. (9) therein) for T,^ can be recast in the form: 

£( 2 ) = C 2 (F)Nd 1 + C 2 (F)n f d 2 + C 2 (F) 2 (d 3 + d i ), (3.4) 
d x = -0.001940(6) d 2 = 0.00237236(16) (3.5) 
d 3 = -0.081429(8) d 4 = 0.01516325(12) (3.6) 

(C 2 (F) is the quadratic Casimir operator in the fundamental representation). In 
the above, contributions proportional to C 2 (F) 2 have been separated into d 3 (arising 
solely from the two diagrams with a tadpole made out of the 4-gluon vertex) and d^ 
(coming from the remaining diagrams). The extension to an arbitrary representation 
R is now immediate: 

£ (2) = C 2 {R)Ndi + 2C 2 (R)T R n f d 2 + C 2 {R)C 2 (F)d 3 + C 2 (R) 2 d 4 , (3.7) 

where the quantities di are left unchanged. 

The prediction from perturbation theory can be improved by resumming a spe- 
cific infinite class of gauge invariant tadpole diagrams. This method is known under 
the name of cactus dressing P3J , and it has been shown to provide improved estimates 



for various quantities of interest, bringing them closer to nonperturbative results. 
Unlike other approaches for improving perturbation theory, such as Refs. j|6], fK] 



this approach does not rely on any Monte Carlo data as input, and it is therefore 
ideally suited for an exploratory study, such as the present one. Cactus resummation 
for the one-loop result of the critical mass simply amounts to dividing the result by 
a factor cq (denoted (1 — w(go)) in Ref. |2l|), which is independent of the fermion 



representation (since it arises from an all-order resummation of gluon diagrams), but 
depends on N and on the bare coupling constant go- The factor Co is the solution of 
the following equation: 



ue -u(N-l)/(2N) 



N-l i / \ n , s 

L 1 N _ 1 (u) + 2L 2 N _ 2 (u) 



9 2 (N 2 -l) _ 9l ( 



N 

[Lp are Laguerre polynomials). For N = 2 and N = 3, Eq. (|3.8j ) simplifies to: 

,2 



c = e 

,2 



-0o/(16co) A_|L^) ) (N = 2) 

co^e flg/(12 ~Co) (l-i- Q + 3^), (iV = 3). 



(3.9) 



Figure ^ presents c |^| as a function of g 2 , for N = 2 and N = 3. The range of 
c/q values, for which a solution exists, extends from g 2 = (where Cq = 1) up to 
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1.10 



1.00 



0.90 



o° 0.80 



0.70 - 



0.60 



0.50. 




0.0 0.2 



Figure 4: cq as a function of g$ for N=2 (red line) and N=3 (green line). 



16/y9e ~ 3.23 (N = 2) and 1.558 (N = 3); this covers the whole region of physical 
interest. The one-loop resummed result thus simply yields: 



m. 



9 2 2C 2 (R) 



(i) , (i) 

1 + c 2 



/co, 



(3.10) 



where [c^ + cij ] is given in Eq. Q3.3| ) and c may be read off Fig. f|. 

In actual numerical simulations, the hopping parameter k is used instead of 
the bare mass itlq. Based on the studies in Ref. ||45|| , we decide to use the one- 
loop resummed result as our estimate of the location of the massless theory in the 
space of bare parameters. The critical value of the hopping parameter is given by: 
k c = l/(2m c + 8), where m c is as in Eq. ( |3.10|) . 

In order to compute the decay constant of pseudoscalar mesons, the renormaliza- 
tion of the axial current Za is needed. While a nonperturbative determination of the 
renormalization constant is desirable, for the first exploratory studies we shall again 
rely on perturbation theory to determine Za- The accuracy of perturbation theory 
can be estimated by comparing the results for QCD, where both a perturbative and 
a nonperturbative computations are available. 

The perturbative renormalization for the axial current with Wilson fermions was 
originally computed in Ref. E3] - see also Ref. [4{| for a useful collection of results 



for the numerical integrals that appear in lattice calculations. The computation 
of the renormalization constant for fermion bilinears requires the computation of 
the vertex functions, and of the fermion wave-function renormalization. For both 
quantities the Feynman diagrams that appear at one loop depend on the quadratic 
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Casimir of the fermion representation. They are therefore readily generalized to an 
arbitrary representation, e.g. 

v - 1 9o C 2 (i?)20.62, (3.11) 



167T 

Z A = l~ ^C 2 (R)15.7, (3.12) 

where the numerical factor is determined by numerical integrals that do not depend 
on the fermion representation. Using the values for C 2 (i?) in the appendix, the values 
of Za and Z v can be easily computed for a generic representation. 

The one- loop perturbative results for Zy, Z^ can also be improved via cactus 
dressing; such an improvement has been known to work rather well with fermions in 
the fundamental representation f50| . To this effect, all that is required is a substitu- 
tion of #o by go/c in Eqs. ( |3.11| , |3.12|) , where c is again read off Fig. f|. 



4. Generalized Aoki phases 

The low-energy dynamics of the pseudo Goldstone particles is determined by the 
pattern of chiral symmetry breaking. For theories with Uf > 1 and fermions in 
a complex representation, such pattern is the same as for QCD, and therefore we 
expect the same Lagrangian to define the dynamics of the effective theory; clearly 
the low-energy constants that appear in the Lagrangian do depend on the specific 
model under study. In this section we discuss the general structure of the chiral 
Lagrangian for fermions in arbitrary representations, and the possible phases of the 
theories discretized on the lattice. 

4.1 SU(2) x SU(2) -► SU(2) 

For the theory Tl that we have been considering in this work, we have two flavours in 
a complex representation of the gauge group and therefore the usual SU(2) x SU(2) 
chiral Lagrangian is expected to determine the dynamics in the low-energy theory. 



Following the notation used in Ref . |5lJ , and including the symmetry breaking terms, 
yields the Lagrangian: 

C = ^Tr (d»Yi%Tl) + jTi (S + £+) - g {Tr (S + E*) } 2 . (4.1) 

The pion decay constant in technicolour theories is related to the vacum expectation 
value of the (composite) Higgs field, which yields F ~ 250 GeV. Lattice artefacts 
for Wilson fermions enter in the coefficients of the symmetry breaking terms: 

d ~ A 3 (m + aA 2 ) , c 2 ~ A 2 (m 2 + maA 2 + a 2 A 4 ) , (4.2) 

where A is again the hadronic scale for the theory under consideration. The pattern 
of symmetry breaking depends on the coefficients C\ and c 2 . Since these coefficients 
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depend on the PCAC mass m and the lattice spacing a, the phase diagram can be 
mapped into the plane of the bare parameters mo, go, used in lattice simulations. 

The analysis in Ref. [SI] remains unchanged; for the theory with two flavours 
the field S can be parametrized as: 

S = A + iB ■ a, A 2 + B 2 = 1, (4.3) 

and the potential becomes: 

- Cl A + c 2 A 2 , (4.4) 

so that the minimum of the potential So = Aq + zBo ■ a can develop a non-trivial 
B only if \Aq\ < 1. For c 2 > a region of width ~ (aA) 3 may exist, where the 
minimum of the potential leads to an Aoki phase. Hence the approach to the chiral 
limit in theories with two Wilson fermions in any complex representation is similar 
to the one observed in QCD: as the quark mass is reduced at fixed lattice spacing, 
flavour symmetry is broken and two massless Goldstone bosons appear. The actual 
values of c\ and c 2 depend on the dynamics of the theory under study, and need 
to be estimated for the cases of interest. Nevertheless, in all cases, the chiral limit 
is entangled with the continuum limit, and the quark mass cannot be lowered to 
arbitrarily small masses at fixed lattice spacing. 

4.2 SU(4) -> SO(4) 

In considering theories in arbitrary representations, different patterns of chiral sym- 
metry breaking may occur. The symmetry breaking patterns and the effective the- 
ories describing the low-energy dynamics in these cases have been studied e.g. in 
Refs. |52|, [53], |54], Using this effective theory framework, we discuss the possibility 



of having an Aoki phase in one case that arises in phenomenologically interesting 
theories, namely the breaking pattern SU(4)— >SO(4). 

This symmetry breaking pattern appears for two Dirac fermions in the adjoint 
representation of the gauge group, see e.g. the theory T2 discussed above. Denoting 
by S the Goldstone matrix, the relevant effective potential for the study of the Aoki 
phase is: 

V = -^Tr [S + St] + ^ {Tr (S + S+) } 2 . (4.5) 

4 lb 

Here S transforms linearly under the global symmetry group SU(4), i.e. 

5 -> g^g T with g G SU(4) , (4.6) 

and 

9 a 

S = S exp(^^X a ). (4.7) 

J 7T 
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In discussing the possibility of an Aoki phase, we are interested in finding at least 
one vacuum configuration where the condensate is not proportional to the identity. 
Indicating with X a and a = 1, ... ,9 the generators spanning the SU(4)/ SO(4) quo- 
tient space [0, Q (an explicit representation of the matrices X a is reported in the 
appendix), we look for solutions of the form: 

9 

S = 2V2(A + tJ2 Aaxa ) > ( 4 - 8 ) 

a 

where A a are real coefficients. With this Ansatz, and using the explicit expressions 
for the generators, the unitarity constraint, i.e. SqSq = 1, implies: 



(Ai + AI + AfXAl + Al + Af 
A\ 



^3 



Substituting the expression for So in the potential we find: 

V = AAl c 2 - 2 Cl A . (4.10) 

For certain values of the potential coefficients C\,C2, there exists a minimum for 
A smaller than one. For instance, by taking A% = A 2 = A 8 = A$ = 0, but a 
nonzero A 3 , we can satisfy the unitarity constraints and SO (4) breaks spontaneously 
to U(l) x U(l). That this is the correct symmetry of the vacuum can be checked by 
determining which SO (4) generator commutes with X 3 . We have checked that the 
S 3 and 5* 4 generators of SO (4) explicitly constructed in Ref. are left unbroken and 
constitute the generators of U(l) x U(l). In this case we would expect the emergence 
of an Aoki phase with four Goldstone bosons associated to the breaking of SO (4) to 
U(l) x U(l). 

4.3 SU(4) -> Sp(4) 

This symmetry breaking pattern does appear for fermions in pseudo-real representa- 
tions. The analysis is similar to the one done before, now using the five (rather than 
nine) X generators presented in the appendix of Ref. We seek for a solution 



using the same Ansatz as in the preceding subsection. The potential evaluated on 
the Ansatz is identical to the one in Eq. ( 4.10| ). In this case the unitarity constraint 



for E yields the condition: 

Al + A\ + A\ + A\ + A\ + A\ = 1 , (4.11) 

For A = 1 one recovers the Sp(4) symmetry. On the other hand a minimum with 
Ao < 1 implies that the Sp(4) symmetry is spontaneously broken to SO (4). We 
hence have four broken generators of Sp(4) corresponding to S l ,S 2 ,S 9 and S 10 in 



the appendix of Ref. [p3|| . Again we find that an Aoki phase is possible, with four 
massless bosons. 
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4.4 Eigenvalues of the Dirac operator 

Finally let us comment briefly on the role of the small eigenvalues of the Dirac opera- 
tor, since they play a crucial role in theories where chiral symmetry is spontaneously 
broken. It was indeed realized long ago that in QCD the chiral condensate is related 
to the density of eigenvalues of the Dirac operator. This property is encoded in the 



Banks-Casher relation 155 



(0\qq\0) = -vrp(O), (4.12) 

where p(X) is the number of eigenvalues in the interval dX per unit volume. 

Following e.g. the derivation in Ref. [pEfl , it is straightforward to show that the 
same relation holds independently of the fermionic representation. As a consequence, 
we expect to have a finite density of eigenvalues around A = for any theory that 
breaks chiral symmetry spontaneously developing a non-zero chiral condensate, and 
hence the number of small eigenvalues of the Dirac operator grows with the four- 
dimensional volume as the mass tends to zero. This phenomenon is directly related to 
the spontaneous breaking of chiral symmetry, and does not depend on the particular 
lattice discretization. 

Moreover, for lattice formulations that do not preserve chiral symmetry, such as 
the Wilson formulation, the spectrum of the hermitian Dirac operator is not bounded 
from below by the bare quark mass. As the quark mass is lowered at fixed lattice 
spacing, the probability of finding an exceptionally small eigenvalue becomes non- 
negligible. These small eigenvalues lead to algorithmic instabilities, violations of 
ergodicity, and sampling inefficiencies, which could seriously distort the output of 
numerical simulations [^]. A careful study of the spectrum of the Dirac operator for 
theories in higher representations is therefore necessary in order to determine a safe 
region for simulating Wilson fermions. This is particularly important as one tries to 
study the phase diagram of novel and unknown theories. 



5. Conclusions 

Gauge theories with fermions in higher-dimensional representations have been put 
forward in several contexts; they are important both for phenomenological and theo- 
retical studies. Some of them provide viable candidates for strong electroweak sym- 
metry breaking, that are not ruled out by precision measurements. On a more the- 
oretical side, they "interpolate" between supersymmetric and non-supersymmetric 
theories, thereby opening new ways to try to tame the nonperturbative dynamics of 
gauge theories. Due to the recent progress in numerical simulations of gauge theories 
on the lattice, it is now possible to simulate these theories, and some first works 
in this direction have already appeared. In this work, we have generalized known 
analytical results to the case of fermions in arbitrary representations. In particular, 
we have considered the scaling as the continuum limit is approached, the location of 
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the critical bare hopping parameter corresponding to massless quarks, the renormal- 
ization of fermionic bilinears, and the phase structure as the quark mass is lowered at 
fixed lattice spacing. The results presented here provide some insight on the unknown 
phase diagram of these theories and will be useful to guide forthcoming simulations. 
Definitive answers on the strong dynamics of such theories, and therefore about their 
viability as phenomenological candidates, can only be provided by actual numerical 
simulations. 
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A. Group— theoretical factors 

The normalization of the generators in a generic representation R of SU(AT) is fixed 
by requiring that: 

[n,T b R ] = if abc T c R , (A.l) 
where the structure constants f abc are the same in all representations. We define: 

tv R (T a T b ) = tr (T*7*) = T R 5 ab , (A.2) 
E FSFrIab = C2(R)6ab, (A.3) 

a 

and hence: 

where dn is the dimension of the representation R. The quadratic Casimir operators 
may be computed from the Young tableaux of the representation of SU(AT) by using 
the formula: 

1 f m 2 \ 

C 2 {R) = -\nN + J2n i (n i + l-2t)- 1 ^\ (A.5) 

where n is the number of boxes in the diagram, i ranges over the rows of the Young 
tableau, m is the number of rows, and rij is the number of boxes in the i-th row. 

The quantities dn, Tr, C^-R) are listed in Table § for the fundamental, adjoint, 
2-index symmetric, and 2-index antisymmetric representations. 
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R 


d R 




C 2 (i?) 


fund 


N 


1/2 


(iV 2 - 1)/(2JV) 


Adj 


N 2 -l 


JV 


JV 


2S 


N(N + l)/2 


(iV + 2)/2 


C 2 (F)2(AT + 2)/(Ar + l) 


2AS 


N(N- l)/2 


(JV-2)/2 





Table 2: Group invariants used in this work 



X 



with 



The generators X a spanning the SU(4)/SO(4) quotient space are defined as: 
1 /V \ 1 / 



2v/2 VO (^f 



1< i < 3; X 1 



D 4 



T 



D 5 — it , D 7 = it 3 



2V2\(Dy 



D 8 = t 1 , 

D 9 = IT 1 



4 < i < 9, (A.6) 



(A.7) 
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